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The formalism of electronic density-functional-theory, with Hubbard-U corrections (DFT+U), is 
employed in a computational study of the energetics of Ui_xCex02 mixtures. The computational 
approach makes use of a procedure which facilitates convergence of the calculations to multiple self- 
consistent DFT+U solutions for a given cation arrangement, corresponding to different charge states 
for the U and Ce ions in several prototypical cation arrangements. Results indicate a significant 
dependence of the structural and energetic properties on the nature of both charge and cation 
ordering. With the effective Hubbard-U parameters that reproduce well the measured oxidation- 
reduction energies for urania and ceria, we find that charge transfer between U(1V) and Ce(lV) ions, 
leading to the formation of U(V) and Ce(lll), gives rise to an increase in the mixing energy in the 
range of 4-14 kj/mol of formula unit, depending on the nature of the cation ordering. The results 
suggest that although charge transfer between uranium and cerium ions is disfavored energetically, 
it is likely to be entropically stabilized at the high temperatures relevant to the processing and 
service of urania-based solid solutions. 



INTRODUCTION 



The thermophysical properties of urania-ceria mix- 
tures have been extensively researched for more than 30 
year^Sm, Interest in this system stems primarily from 
two issues related to the performance of oxide nuclear 
fuel. First, cerium is a fission product in urania nuclear 
fuels that is highly soluble in the fluorite UO2 structur^ 
its effects on phase stability, thermal conductivity and 
ionic diffusion are thus important issues relevant to nu- 
clear fuel performance, particularly for future potential 
applications involving high burnup. Second, urania-ceria 
has been investigated as a surrogate system for urania- 
plutonia mixed-oxide (MOX) fuels, owing to the simi- 
lar oxidation states and ionic radii for Ce and Pu ions. 
In spite of the interest in this system, and the exten- 
sive experimental research that has been devoted to it to 
date, the thermodynamic properties and phase diagram 
of urania-ceria remain incompletely characterized. 

Due to the importance of the solution-thermodynamic 
properties and phase equilibria on the performance of 
nuclear fuels, renewed efforts have been undertaken to 
develop accurate thermochemical models for these ma- 
terials. A prominent approach is the application of 
the so-called CALPHAD (Calculation of Alloy Phase 
Diagrams) methodologj-^, involving the parametrization 
of compound-energy- form alisrrP free-energy models for 
UO2 and its mixtureJsEl, Por these efforts, a challeng- 



ing issue is the lack of complete databases of experimen- 
tal thermodynamic data, as well as the often conflict- 
ing nature of the results that are available. In CAL- 
PHAD modeling of metal alloy systems, a strategy that 
has been effectively pursued to improve the accuracy of 
the thermodynamical models involves the application of 
first-principles calculations as a framework to augment 
available experimental data in the development of robust 
thermodynamic databased^"^. Over the last decade ad- 
vances in first-principles methods have led to widespread 
use of these techniques in the modeling of defect struc- 
tures and energetic properties in pure UOjP^'^. To date, 
however, they have not been applied extensively in the 
study of mixed oxide thermodynamic properties. 

The present work involves application of first- 
principles methods, based on the Density Functional 
Theory with Hubbard-U (DFT-I-U) formalism (see sec- 
tion HA), to the study of mixing energetics in urania- 



ceria solid solutions with stoichiometric oxygen composi- 
tions. The work represents an initial step towards more 
complete modeling of solution thermodynamics for both 
stoichiometric and non-stoichiometric oxygen composi- 
tions in this system. It also provides insight into an is- 
sue that has been extensively discussed in the context 
of the properties of urania-ceria solid solutions, namely 
the presence of mixed charge states for the uranium and 
cerium cations. Specifically, magnetic susceptibilitjff^ 
and X-Ray Photoemission Spectroscopy (XPS) mea- 
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surements have shown evidence of charge-transfer pro- 
cesses by which a fraction of the U(IV) and Ce(IV) 
ions are oxidized and reduced, respectively, to U(V) 
and Ce(IIl|^. However, conflicting conclusions have 
been drawn from X-Ray Adsorption Near-Edge Spectra 
(XANES) measurement where it was concluded that 
Ce and U retain oxidation states of IV in Ce-rich solid 
solutions. The current work presents a computational 
methodology that can be used to systematically investi- 
gate the relative energetics of solid solutions with ideal 
versus mixed-charge states, and suggests a picture in 
which Ui_xCex02 compounds with mixed charge states 
are energetically disfavored while being entropically sta- 
bilized at high temperatures. 

The remainder of this paper is organized as follows: 
we begin by describing the computational methodology 
and follow with an explanation of how we converge to 
and characterize different charge states in the DFT-I-U 
calculations. The following section presents results for 
energetics, ionic structure, and electronic structure, aug- 
mented with classical-pair-potential studies extending 
the work to disordered systems. We conclude by dis- 
cussing implications of the first-principles results for the 
mixing thermodynamics in this system. 



II. METHODS 

To investigate the energetics of charge and cation or- 
dering in the urania-ceria system, we consider a set of 
fluorite-based superstructures in which the cations are ar- 
ranged in accordance with several fcc-based prototypical 
ordered compounds. Cation arrangements in these struc- 
tures are built from ordering waves along the <100>, 
<210> and <111> "special-point" ordering directionJiS, 
and include the so-called "Lifshitz structures" for order- 
ing on the fee cation sublattice of the fluorite structure. 
All structures considered in this work are stoichiomet- 
ric, i.e., the oxygen fluorite sublattice is fully occupied 
with no interstitial oxygen atoms present. Experimen- 
tally, urania-ceria solid solutions are observed to be con- 
figurationally disordered on the cation sublattice. How- 
ever, in the present work we consider only the ordered 
structures shown in Fig. [T] and listed in Table |T] in the 
DFT-I-U calculations. To model the energetics of disor- 
dered solid solutions, we employ pair-potential models 
tuned to the DFT-I-U results, in supercell calculations of 
random mixtures. 

In our calculations we consider two types of ionic 
charge states. In the first, which we will refer to as "ideal 
charge states" (ICS), all cerium and uranium ions possess 
the 4-1- charge state (see discussion below for the defini- 
tion of the charge states based on the results of the elec- 
tronic structure calculations). In the second, which we 
will refer to as "mixed charge states" (MCS), a fraction 
of the cerium and uranium ions are respectively reduced 
and oxidized to 3-1- and 5-1- charge states in a 1:1 ratio of 
Ce(III):U(V) to preserve overall charge neutrality. 



Composition 


Prototype 


Strukturbericht 


Ordering Wave 


AgB 


AuCu3 


LI2 


<100> 


A2B 


MoSi2 


Glib 




A2B2 


AuCu 


Llo 




A2B 


Cdl2 


C6 


<111> 


AB 


CuPt 


Lll 




A3B 


AlgTi 


DO22 


<210> 


A2B 


MoPt2 






A4B4 


NbP 







TABLE I. Compositions, prototypes, Strukterbericht and as- 
sociated ordering wave families for the prototype ordered 
structures considered in the DFT-fU calculations. 



AuCug MoSij AuCu 



<100> 




Cdl2 CuPt 




FIG. 1. (color online) An illustration of the prototype ordered 
structures considered in the DFT-f U calculations for the five 
compositions of mixed oxide explored. Ordering waves are 
identified, and the stacking planes are shown in gray. The 
oxygen cubic sublattice is also shown. 



A. First Principles Calculations 

All DFT-I-U calculations have been performed employ- 
ing the formalism of Dudarev et al^^, as implemented in 
the Vienna ab-initio simulation package (VASP)^*' In 
the implementation used in this work, we have empl oyed 
the projector augmented-wave (PAW) methocP^EH and 
the Perdew-Becke-Ernzerhof generalized gradient ap- 
proximation (PBE-GGA)25i26] The PAW-PBE poten- 
tials employed in this work are those designated "U", 
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"Ce" , and "O" in the VASP library. The electronic wave 
functions were expanded in a plane-wave basis set with a 
cutoff of 500 eV, and the electronic states were sampled 
using k-point meshes centered on the origin, with a den- 
sity equivalent to that of a 4x4x4 mesh in the Brillouin 
zone of the ideal fluorite structure. The structures were 
fully relaxed with no symmetry constraints imposed, un- 
til the magnitude of the forces was below 10"^ eV/A and 
the pressure was below 1 kBar. With these settings the 
calculated mixing energies are estimated to be converged 
to within 1 meV per (U,Ce)02 formula unit (FU). The 
mixing enthalpy is defined as follows: 

Hmix — Hu^Cei-^02 ~ {xHu02 + (1 ^ x)Hce02)i (1) 

where Hij^cei^^02 denotes the energy (per cation) of a 
urania-ceria mixture with x denoting the cation mole 
fraction of uranium ions, and similarly for Huo2 ^^'^ 
Hce02 ■ -'^U of the results presented in this manuscript 
were calculated using the scalar-relativistic approxima- 
tion. For the LIq ICS and MCS structures, we also 
performed calculations including spin-orbit coupling and 
the resulting mixing energies were found to change by a 
magnitude of less than 3.5 meV/FU for the ICS, and 25 
meV/FU for the MCS. 

The differences in energy between MCS and ICS for a 
given structure are highly sensitive to the choice of the 
parameter U^g-^U-J in the DFT-I-U formalism employed 
in this work. The values of U^g- for U 5/ and Ce 4/ 
electrons utilized in the present calculations yield a close 
match between calculated and experimental values of rel- 
evant oxidation and reduction energies. For uranium, the 
choice of Ugg:=3.99 eV, which was proposed originally 
by Dudarev et al^ and has been widely used in the lit- 
erature, was found in the present study to give rise to 
good agreement between calculated and measured values 
of the enthalpy changes for the following two oxidation 
reactions: (i) UO2 + 5 O2 ^ 7-UO3, (ii) 3 UO2 + O2 
— ?► a-U30g. Descriptions of 7-UO3 and the a-U30g can 
be found in Refs.^ ancPll^ respectively. With \]qQ—3.99 
eV, and employing the correction for over-binding of the 
oxygen molecule discussed in Ref.'^, we obtain values 
of A IIrxn=-l-45 eV and -3.44 eV for reactions (i) and 
(ii) , which agrees well with the corresponding experimen- 
tal values of -1.44 eV and -3.31 eV, as determined from 
the CODATA Key Values for Thermodynamic^. For 
cerium we choose the value Uggf=3 eV, which was shown 
by Andersson et al^ to yield good agreement between 
calculated and experimental values for the enthalpy of 
the reduction reaction: Ce203 -f 5 O2 — >■ 2 Ce02. 



B. Electronic and Ionic Relaxation 

An important issue that has been recently discussed in 
the literaturcPtSll concerns the propensity for DFT+U 
calculations to converge to multiple self-consistent solu- 
tions corresponding to different orbital occupations. To 



ensure that such calculations converge to solutions that 
are the lowest-energy electronic states (or near them), 
several different methods have been proposed. In the cur- 
rent work we employ the approach described by Meredig 
et al^, which involves a slow localization of the / elec- 
trons in a series of DFT-f U calculations with incremen- 
tally increasing values of U^g. Specifically, one begins by 
performing a GGA (with Ugff=0) calculation until ionic 
positions and charge density are converged. The wave- 
function eigenvalues, charge densities, and atomic posi- 
tions are then used as the starting point for a calculation 
with a small non-zero value of Ugg-, again performed to 
convergence of ionic position and charge density. The re- 
sults of this calculation are then used as the starting point 
for a calculation with an incrementally higher Ug-g:, and 
this process is continued until the desired Ug^ is reached. 
Steps in Ugg- of approximately 0.1 eV have been found 
sufficient to reproducibly converge the low energy struc- 
ture of UO2 • 

For urania-ceria mixtures, this method for converg- 
ing the electronic structure becomes somewhat ambigu- 
ous. Specifically, Ugg can be incremented for each ion 
type simultaneously or in two separate stages. However, 
this ambiguity has been exploited in the current work 
to enable convergence of the electronic structure to the 
ICS and MCS charge configurations. Indeed, we found 
that the ICS charge state (i.e., composed of Ce(IV) and 
U(IV)) is obtained if Ug^ is first ramped up to its final 
value on uranium 5/ electrons, while holding the value of 
UgQ-^O for the cerium 4/ electrons. Alternatively, if Ugg- 
is first applied to cerium 4/ electrons, the MCS charge 
state (i.e., containing Ce(III) and U(V)) is obtained. 

The oxidation states of the ions resulting from the cal- 
culations described above are identified based on three 
considerations: the interatomic bond lengths, the num- 
ber of states in the occupied electronic partial densities of 
states, and the local magnetic moments. For structures 
with Ce(III) ions, the number of states in the occupied 
Ce 4/ projected densities of states is found to integrate 
to approximately one electron (cf. Fig. [4]), while this 
band lies above the Fermi level for structures containing 
only Ce(IV) ions. Similarly, the number of states in the 
occupied U 5/ projected densities of states is approxi- 
mately two and one for U(IV) and U(V), respectively. 
Corresponding to these different charge states are dis- 
tinct values of the local magnetic moments. For U(IV), 
U(V), Ce(III) and Ce(IV) ions the calculated local mo- 
ments are close to the ideal values of 2, 1, 1 and 0, re- 
spectively, as expected from Hund's rule, which assumes 
that these moments originate from the single and paired 
/ electrons. A final consistency check associated with the 
labeling of the charge states is made based on the calcu- 
lated relaxed bond lengths. Oxidation of U(IV) to U(V) 
is found to lead to the expected decrease in the average 
relaxed U-0 bond lengths, while reduction of Ce(IV) to 
Ce(III) leads to an increase of the Ce-0 bond lengths (cf. 
Fig. ^. 
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C. Classical Pair Potentials 

To understand the implications of the first-principles 
results for the thermodynamic properties of disordered 
solid solutions, we have used the energetic results pre- 



sented in section III to optimize ionic pair-potential mod- 
els for urania-ceria solid solutions. Specifically, we em- 
ploy the polarizable shell-model potentials of the Buck- 
ingham type for O, U(IV), Ce(III) and Ce(IV) published 
in RefspSffill^ ^ similar potential for U(V) has been previ- 
ously developed by one of us'^^ . The interaction parame- 
ters for this potential are A=2386.42 eV, p=0.3411 A and 
C=0 eV-A^. For U(V), an important issue concerns the 
magnitude of the fifth ionization potential, which cannot 
be calibrated to any direct experimentally measurable 
properties, but is required in order to compare energetic 
differences between ICS and MCS states. Calculations 
by Pyper and Grant estimate this value at 46.57 e\^. 
We found that an optimal agreement between potential- 
model predictions and the DFT-I-U energetics is obtained 
with a value of 47.62 eV, an increase of approximately 
2%. Using this value for the fifth-ionization energy, and 
the published parametrizations for the other cations, the 
classical-potential models match the DFT-I-U mixing en- 
ergies of X = 0.25, 0.50, and 0.75 in Fig. [2] (discussed in 
more detail in the Results section) to within a standard 
deviation of 4.5 kJ/mol. 

With the pair potentials, the focus on urania rich dis- 
ordered solid solutions is motivated, in part, by ongoing 
experimental work in this systeni-^. It is also within this 
regime that the solution thermodynamics of such materi- 
als is relevant to nuclear fuels in particular. The calcula- 
tions employing ionic pair potentials proceeds as follows. 
Supercells containing 256 cations were constructed with 
random occupations of the cation sublattice at composi- 
tions of 25, 21.9, 15.6, 12.5, 6.25, and 3.1 cation percent 
Ce. For each of these cation configurations, we also con- 
sidered ICS states as well as two MCS states with 50 
and 100 % of the Ce ions reduced to the a charge state 
of III (and a corresponding number of U ions oxidized 
to a charge state of V to maintain charge neutrality). 
At each of the resultant 18 compositions, 300 supercells 
were generated with different random-number seeds for 
each composition. For each of the constructed supercells, 
full geometry relaxations were performed employing the 
GULP softwarcP. 



III. RESULTS 

A. First-Principles Results 

Figure [2] plots DFT+U calculated values of the mix- 
ing energies for each of the prototype ordered structures 
shown in Fig. [T] for both ICS and MCS states. Con- 
sidering first the ICS structures, the mixing energies are 
relatively small in magnitude, ranging from very weakly 



negative values of 1 kJ/mol to positive values less than 
5 kJ/mol in magnitude. The MCS structures are cal- 
culated to be higher in energy than their corresponding 
ICS analogues at the same composition. However, the 
difference in energy between ICS and MCS states is seen 
to be sensitive to the nature of the cation ordering. For 
example, focusing on the structures at equiatomic com- 
position, the AuCu and NbP structures, which are built 
from <100> and <1^0> special-point ordering waves, 
respectively, have energies that increase by 14 kJ/mol 
in going from the ICS to MCS states. By contrast, the 
CuPt structure, built from < 5 55 > ordering waves, is 
only 4 kJ/mol higher in energy in the ICS relative to the 
corresponding MCS state. A comparison of the results 
at each composition also shows an interesting trend that 
the structures with the highest energies in the ICS states 
have the lowest energies in the MCS states. In general, 
the variation in the energy differences between ICS and 
MCS states for the different structures suggests a sig- 
nificant coupling between the nature of the cation order 
and the energetics associated with charge transfer in the 
system. 
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FIG. 2. Zero temperature mixing enthalpies for ordered struc- 
tures as a function of Ce fraction on the cation sites. Energies 
for MCS compounds are shown with gray squares, while those 
for ICS compounds are indicated by black circles. 

In Fig. [3] we compile results for the calculated atomic 
structures of the ICS and MCS compounds. The re- 
sults are plotted in the form of distributions of bond 
lengths between the different cation species and their 
eight nearest-neighbor oxygen ions. These distributions 
are derived from the fully relaxed structures of all the 
ICS and MCS compounds considered in the DFT-I-U cal- 
culations. For ICS ordering, the averaged lengths for 
U(IV)-0 and Ce(IV)-0 bonds are 2.393 and 2.385 A, re- 
spectively, reflecting a slightly smaller size for the Ce(IV) 
ions relative to U(IV) in the fluorite structure. In the 
MCS compounds, the U(IV) and U(V) bond-length dis- 
tributions overlap significantly, although the former show 
a clearly smaller average bond length, as expected: the 
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FIG. 3. Bond length distributions for ICS ordering (left panel) 
and MCS ordering (right panel). Distributions are summed 
over all structures considered in the DFT+U calculations, 
amounting to a total of 384 metal-oxygen bonds for the ICS, 
and 352 for the MCS compounds. 



averaged values for the U(IV)-0 bond length in the MCS 
compounds is 2.399 A, which is similar to its averaged 
value in the ICS structures, while the corresponding value 
for U(V)-0 in the MCS compounds is 2.336 A. The dis- 
tribution of Ce-0 bond lengths in the MCS compounds 
is seen to be significantly narrower than that of the U- 
O bonds in the same structures. The average Ce(IV)- 
O bond length in the MCS compounds has a value of 
2.386 A which is very similar to that in the ICS struc- 
tures. The Ce(III)-0 bond length is larger, as expected, 
with an average value in the MCS compounds of 2.458 
A. Overall, the results display the expected trends as- 
sociated with changes in the cation radii with different 
oxidation states. The broader distributions for the U-0 
bond lengths in the MCS relative to the ICS compounds 
suggests a higher degree of strain energy in the former 
structures, with the preferred bond lengths for the differ- 
ent cation charge states being accommodated to differing 
degrees depending on the nature of the cation ordering. 

The differences in the electronic structures for ICS ver- 
sus MCS compounds is illustrated by the calculated elec- 
tronic densities of states (DOS) for the equiatomic AuCu 
prototype structure in Fig. |4] In this figure, the ICS 
structure features a narrow band just below the Fermi 
level corresponding to the occupied 5/ orbitals for the 
U(IV) ions. This band is split off from the oxygen 2p 
band with a gap of 0.80 eV. In the MCS compound, the 
U 5/ states shift to significantly lower energies with the 
highest density of states for these orbitals found at the 
edge of the O 2p band. The U(V)-0 bonding thus dis- 
plays a significantly stronger degree of covalent character, 
characterized by a stronger hybridization between the U 
5/ and O 2p states in the MCS compound. The occupied 
Ce 4/ states for the Ce(III) ion in the MCS compound 
are observed to be split off from the U 5/ and O 2p bands 
by a gap of approximately 2 eV, and these states display 
a much lower degree of hybridization with the O 2p or- 
bitals. 




-6 -4 -2 2 4 -6 -4 -2 2 




FIG. 4. (color online) Total density of states for AuCu or- 
dered: a) ICS; and b) MCS. Partial density of states for AuCu 
ordered: c) ICS; and d) MCS. 



B. Pair-Potential Modeling 

The DFT-I-U results presented in Figs. [2][6] have been 
derived for long-range-ordered compounds featuring pro- 
totypical ordering of U and Ce ions on the fee cation 
sublattice. As discussed above, stoichiometric urania- 
ceria mixtures observed experimentally are disordered 
solid solutions, lacking measurable cation order. To 
investigate the implications of the DFT-fU results for 
the mixing thermodynamics of disordered solid solutions 
we have undertaken calculations based on classical ionic 
pair-potential models, as described in section [ll C| 

Figure [5] plots the averaged values (symbols) and stan- 
dard deviations (indicated by error bars) of the mixing 
energies calculated by averaging over supercell configu- 
rations with randomly generated cation and charge dis- 
order. The mixing energies for the ICS states are small 
and positive in magnitude, consistent with the DFT-fU 
results for the ordered ICS structures in Fig. [2j With 
increasing degree of charge transfer, the mixing energy 
is seen to increase and is analogous to the DFT calcula- 
tions, in which charge transfer is disfavored energetically. 
Finally, the increased variation in the energies of struc- 
tures with high degrees of charge transfer is reflected by 
the larger standard deviations in the calculated mixing 
energies. 

Although charge transfer is found to lead to a higher 
energy of the MCS relative to the ICS states in Figs. [2] 
and[5j the energy difference is relatively small, such that 
the MCS configurations are expected to be sampled for 
entropic reasons at the temperatures used in the syn- 
thesis and applications of urania-based fuels. To demon- 
strate this point, we show in Fig. [6]the mixing free energy 
derived at T=1500 K, obtained by combining the cal- 
culated mixing enthalpies with an ideal mixing entropy. 
Specifically, the free energy plotted in the left panel is 
obtained by minimizing the expression: 



6 



go 4 



OC= 1.0 




■ C = 0.5 






AC = 0.0 

































0.03 0.05 0.07 0.09 0.11 0.13 0.15 0.17 0.19 0.21 0.23 0.25 
Cation Fraction of Cerium 

FIG. 5. Average mixing enthalpies obtained from ionic pair- 
potential calculations for structures enumerated by randomly 
populated supercells. The average mixing energy obtained for 
each composition is given by a filled symbol, and the solid lines 
are polynomial fits, included to highlight the trends. Error 
bars represent ±1 standard deviation in all of the data for a 
given composition. "C" is the mole fraction of Ce which is in 
the Ce(III) state. 



A (Xceilll) , XCe(IV) , ^UilV) , Xu(V) ) 

+ RT (2) 

_ i 

with the mass and charge neutrality constrains: 



IV. SUMMARY AND DISCUSSION 



In this work, we have examined the mixing energetics 
of Ui_xCex02 using DFT-I-U and classical pair poten- 
tial methods. Results obtained for prototypical ordered 
structures on the cation sites result in largely structure 
independent, nearly ideal mixing enthalpies for systems 
in which U and Ce ions maintain the ideal 4-1- charge 
state. With the values of Ugg- used in this work, DFT-f U 
calculations show that charge transfer between uranium 
and cerium sites (resulting in U(V) and Ce(III) cations) 
leads to an increase in mixing energy. We note that the 
sign and magnitude of the effect of charge transfer on 
the mixing energetics is sensitive to the choice of the 
Uggf parameters in the DFT-I-U calculations. The re- 
sults obtained in this work are based on values that have 
been found to lead to a good match (accounting for over- 
binding of the oxygen molecule) between calculated and 
measured values for relevant oxidation and reduction re- 
actions in urania and ceria. In the present calculations, 
the magnitude of the increase in AH^j^j^j. resulting from 
the formation of mixed-charge states is highly sensitive 
to the ordering of the cations. Despite this tendency 
for charge transfer to destabilize the mixtures energeti- 
cally, the energy difference between ICS and MCS states 
is small enough that charge disordering may be entrop- 
ically stabilized at high temperatures (1500 K) relevant 
to synthesis and applications of urania solid solutions. 
This point was illustrated in section [ill B| with the aid of 
pair-potential models. 



^ Xi = XCe(IH) + XCe(IV) + Xu(IV) + Xu(V) = 1 (3) 



^xceiiii) + 4:Xce{iv) + 4a;t/(/y) + 5xu(v) = 4 (4) 

where Xj is the fractional cation lattice occupation of 
species i. The free energy surface thus obtained is plotted 
in the middle panel of Fig. [6| 

The degree of charge-transfer is defined as C = 
xce{iH)/{xce{iii) + Xce{iv)), and its value which mini- 
mizes AGjj-^j-jj is shown in the right panel of Fig. |6] The 
corresponding value of the mixing free energy is shown in 
the left panel of the same figure. The results display an 
appreciable degree of charge transfer in dilute solutions 
that decreases with increasing Ce concentrations. This 
behavior reflects the fact the degree of charge transfer is 
driven by entropic contributions that increase with in- 
creasing electronic disorder. We note here that the en- 
tropic stabilization of mixed-charge states, suggested in 
the present calculations for urania-ceria solid solutions, 
has also been shown in recent calculations for LiFeP04 
battery material^^. 



The small differences in energy between the ICS and 
MCS states obtained in this work suggest that the de- 
gree of charge transfer in urania-ceria solutions is likely 
to depend sensitively on oxygen stoichiometry. Indeed, 
it has been shown in XPS studies-'-^ that the concentra- 
tion of Ce(III), U(V), and U(VI) are all strongly depen- 
dent on oxygen content. Notably, however, it is also 
possible that the degree of charge transfer could be dif- 
ferent in the near-surface regions probed by XPS. This 
would provide a possible reason for why charge trans- 
fer is found in XPS studies, while no evidence for mixed 
charge states was obtained by XANES measurement!^. 
At present the most compelling evidence for charge trans- 
fer in bulk samples is that derived from magnetic suscep- 
tibility measurements^^. In light of the present compu- 
tational results, and the existing discrepancies in conclu- 
sions drawn from different experimental methods, mul- 
tiple experimental methods applied to the same set of 
samples would be highly desirable. Given that the cation 
interactions and oxygen stoichiometry are likely to be 
strongly coupled in this system, attention to characteriz- 
ing composition and homogeneity would be particularly 
important for such future studies. 
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FIG. 6. (color online)The center graph depicts the ideal mixing free energy as a function of the two composition variables, 
cerium cation fraction ("x") and the percentage of the Ce ions which are in a charge state of III ("C"). The left- and right-most 
panels give the values of the mixing free energy and the values of the charge transfer parameter ("C"), respectively, which 
correspond to the convex hull of this plane. These illustrate the lowest energy states at each total cerium composition. 
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